Ktz form factors, final state interactions and K ti^ti^ 

decays' 

D. R. Boito and R. Escribano 

Grup de Fisica Tedrica and IFAE, Universitat Autdnoma de Barcelona, 
E-08193 Bellaterra (Barcelona), Spain 

Abstract. We present a model for the decay £)+ — > Tt^Tt^ . The weak interaction part of this reaction is described using 
the effective weak Hamiltonian in the factorisation approach. Hadronic final state interactions are taken into account through 
the Kn scalar and vector form factors fulfilling analyticity, unitarity and chiral symmetry constraints. Allowing for a global 
phase difference between the 5 and P waves of —65°, the Dalitz plot of the £)+ K^k^k^ decay, the Kti invariant mass 
spectra and the total branching ratio due to 5-wave interactions are well reproduced. 
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INTRODUCTION 

In 2002, the analysis of — > K ;r+ ;r+ decays performed by the E79 1 collaboration revealed that approximately 50% 
of these decays proceed through a low-mass scalar resonance with isospin 1/2: the /Tq (800), also called the k[1]. More 
recently, the K^n^n^ decay was revisited by E791 [2] and two other experiments produced analyses based on 

larger data samples, namely FOCUS [3, 4] and CLEO [5]. The main conclusions of the pioneering E791 work have 
been confirmed in both cases. 

In the past, many analyses of Kn scattering data had already claimed the presence of the K pole in the scattering 
amplitude [6, 7, 8, 9]. The most precise and model independent determination of its position in the second Riemann 
sheet was produced in Ref. [10]. Using Roy's equations for Kn scattering [11] and Chiral Perturbation Theory 
(ChPT) [12] Descotes-Genon and Moussallam found 658 ± 13 MeV and = 557 ± 24 MeV [10]. 

Although the experimental results are sound and the K pole is at present theoretically well known, a comprehensive 
and successful description of the reaction D+ — >^ K^n'^n'^ is still not available (for a recent review see Ref. [13]). 
Experimentalists, for the want of a better framework, commonly fit their data with the isobar model which consists 
of a weighted sum of Breit-Wigner-like propagators. Often, a complex constant is added to the amplitude in order to 
account for the non-resonant decays. It is known, nevertheless, that the adoption of Breit-Wigner functions to describe 
the effect of scalar resonances is problematic. 

In the present work, we follow the general scheme where a factorised weak decay amplitude is dressed with FSIs 
by means of non-perturbative Kn form factors. For the weak vertex, we employ the effective weak Hamiltonian of 
Refs. [14, 15] within naive factorisation. The weak amplitude thus obtained receives contributions from colour-allowed 
and colour-suppressed topologies. In the latter, the Kn form factors appear manifestly and the construction of the final 
state is straightforward. The colour-allowed topology is more involved but, assuming the decay to be mediated by 
resonances as suggested by the experimental results, the FSIs in this case can also be written in terms of Kn form 
factors [16, 17]. Therefore, in our description the hadronic FSIs are fully taken into account by the Kn scalar and 
vector form factors. Both form factors have received attention in recent years and are now well known in the energy 
regime relevant to — > K^n^n^ decays. The scalar component was studied in a framework that incorporates all 
the known theoretical constraints in Refs. [18, 19, 20]. The results were subsequently updated and we employ in this 
work the state-of-the-art version given in Ref. [21]. The vector form factor, in its turn, can be studied in KnVx 
decays [22, 23, 24, 25], where the kinematical range is very similar to the one considered in this paper Here we employ 
a description which fulfils analyticity constraints and that was successfully fitted to the Belle spectrum in Ref. [25]. 

This contribution is based on a recent paper where all the issues discussed here are presented in more detail [26]. 
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THEORETICAL FRAMEWORK 



Our phenomenological description of the process £)+ A'^7r+7r+ is based on the effective weak Hamiltonian. At 
the quark level, the decay — >• K~7t'^7t'^ is driven by the transition c sud, i.e. four different quark flavours are 
involved. In this case, only the two tree operators of the weak Hamiltonian have to be taken into account. The amplitude 
forZ)+ — >^ A'~7r+7r+ is given by the matrix element {K^n^n^\J^s\D'^). We assume the factorisation approach to hold 
at leading order (in Aqcd/wic and a,) and as a consequence the amplitude is written in terms of colour allowed and 
suppressed contributions, ^ and ^ respectively, as 

Gf 

si/{D+ ^K~n+n+) = ^cos^ 0c(aiM +^2^) + (?rf ^ tTj^) 

= ^cos^dc[ai{K-K+\st{l-Y5)c\D+){T4\uy^{\-y5)d\Q) 

+a2{K-K+\s'f{^ - Y5)d\0){n+\uY^,{\ - Y5)c\D+)] + {k+ ^ k+) , 

where the last term accounts for the presence of two identical pions in the final state. For the QCD factors 012 we 
use the phenomenological values ai = 1 .2 ± 0. 1 and = —0.5 ± 0. 1, obtained from different analyses of two-body D 
meson decays [27]. The non-perturbative hadronic matrix elements in Eq. (1) involve several Lorentz invariant form 
factors [26]. The ampUtude ^ reads 
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where the Mandelstam variables are defined as m^Km = (pk + Piri)^, "^^ttj ^ ^P'^ ^ P^^i)'^ ^^'^ ^\in2 = (P^i + P^2) 
with + ni\j^ + rr^^^ = mj) + mj. + 2m^. Here, we follow Ref. [16] and write the colour aUowed amplitude 



in terms of the scalar and vector Kn form factors as 

,.cff /.„2 2 \ cA'ff 



where Xs^ and two free parameters that are fixed from experimental branching ratios. 



NUIMERICAL RESULTS 

Let us now investigate in detail the numerical results for our final model which includes the contribution of both ^1 
and ^2 topologies. The corresponding expressions are given in Eqs. (2) and (3). We begin by considering the 5-wave 
description which is, in our opinion, the main aspect of the problem. In our model, the 5-wave FSIs are described by 
the Kn scalar form factor of Ref. [21] in a quasi two-body approach, i.e., we assume that the Kn pairs in Eq. (1) form 
an isolated system and do not interact with the bachelor pion. Moreover, the form factor of Ref. [21] is obtained from 
dispersion relations that fix its phase to be the scattering one within the elastic region [18]. Consequently, our 5- wave 
amplitude has the Kn 7=1/2 scattering phase up to roughly 1.45 GeV where the Krj' channel starts playing a role. 
We compare in Fig. la the experimental results from Refs. [2, 5, 4] with the phase of our S wave. Since we are dealing 
with a production experiment, a global phase difference is expected as compared with scattering results [3]. Therefore, 
we allow for a global phase shift a in our 5-wave amplitude x/siff-^m ^^kji2^ ~^ •'^('"a:;ii '"^A:;r2)- ^^?>- ^® 
also plot as the dot-dashed line the phase of our amplitude shifted by a = —65°. With this shift, we see that up to 1.5 
GeV CLEO's results and ours share a remarkably similar dependence on energy. Inspired by the inspection of Fig. la, 
we consider as our final model the one given by Eqs. (2) and (3) with a shift of a = —65° in the 5-wave phase. 

In order to compare the absolute value of our S wave amplitude with experimental data, we need fix the only two 
free parameters that occur in our model, namely the normalisation constants xf^ and Xv^ ■ The constant Xv^ is fixed in 
order to reproduce the value of the sum of all vector submodes. Then, we fix the scalar normalisation xf^ requiring the 
total branching ratio from our model to match the world average. Taking the central values for a\ and ai this procedure 
gives xf^ = 4.9 ± 0.4 GeV~^ and Xv^ = 4.4 ± 0.6 GeV~^ We can now compare the absolute value of our 5-wave 




FIGURE 1. (colour online), (a) 5-wave phases from the analyses of Refs. [2, 5, 4]. The solid line is the phase of our S-wave 
amplitude with a = 0°, whereas the dot-dashed line is the S-wave phase with a = —65°. The dashed line delimits the Krj' threshold, 
(b) Absolute value of the 5 wave measured in Refs. [2, 5, 4] compared with our model. The amplitudes are normalised according 
to the text. 



amplitude with experimental results. However, since in isobar-like analyses the fit is sensitive only to the relative 
weights of the amplitudes, in order to compare the measurements with our result we need perform a normalisation. 

1 /2 

We define a normalised 5-wave amplitude by S!/g""" {mj.^^ ,mj^j^^) — s^^l /® '^^k%x'^^kti2 I"*^!^) ■ """^^^ amplitude, 
by construction, is free of any global constants that appear in s^^ and has dimension of [Energy] Interpolating 
the results from the tables found in Refs. [2, 5, 4] we can calculate the normalised S wave for each experiment. We 
repeated the same procedure for our total 5-wave amplitude. The results for the S wave are compared with our model 
in Fig. lb. 

We can also perform a Monte Carlo (MC) simulation to obtain a Dalitz plot from our model and compare the 
diagram and its projections with experimental results. For the lack of a true data set, we resort to a MC simulation of 
the original E791data [1]. The obtained diagram is shown in Fig. 2a. Then we performed the same exercise for our 
model and the result is shown in Fig. 2b. Finally, in Fig. 3 we show the projections of the diagrams of Figs. 2a and 2b. 



SUMMARY AND DISCUSSION 

We have presented a model aimed at describing the decay £)+ — > K^%^%^ . The weak amplitude is described within 
the effective Hamiltonian framework with the hypothesis of factorisation. The K% hadronic FSIs are treated in a quasi 
two-body approach by means of the well defined scalar and vector K% form factors, thereby imposing analyticity, 
unitarity and chiral symmetry constraints. We used the experimental values for the total and P-wave branching ratios 
to fix the two free parameters in the model. The relative global phase difference between the S and P waves was fixed 
phenomenologically using the experimental results of Ref. [5]. 

The use of the K% scalar form factor is shown to provide a good description of the 5-wave FSIs. Both the modulus 
and the phase of our S wave compare well with experimental data up to < 1-5 GeV. It is worth mentioning that 
the form factor we used has a pole that can be identified with the k. Furthermore, the model is able to reproduce the 
experimental fit fractions and the total 5-wave branching ratio. Finally, the Dalitz plot arising from the model agrees 
with a MC simulated data set. 

The main hypotheses of our model are the factorisation of the weak decay amplitude and the quasi two-body nature 
of the FSIs. Therefore, the success of our description for niKK < 1 -5 GeV suggests that, in this domain, the physics of 
the decay is dominated by two-body Kn interactions. We are led to conclude that effects not included in our model 
such as the 7 = 3/2 non-resonant K% S wave, the non-resonant 1 — 2 7i;+7r+ interactions and genuine three-body 
interactions, could be considered as corrections to the general picture described here. 

Part of the discrepancy observed in our Dalitz plot is due to the disaccord of our 5-wave amplitude for ihkk > 
1.5 GeV. A possible cause for this disagreement is the fact that factorisation in a three-body decay is expected to break 




FIGURE 2. (a) Monte Carlo simulation for the Dalitz plot of the E791 original analysis [1] (b) Same for our model with a global 
shift of —65° degrees in the 5-wave phase (see text and Fig. 1). The number of independent events is 14185, which correspond to 
the estimate of the signal events in Ref. [1]. 
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FIGURE 3. (colour online). Projections from the MC generated Dalitz plots of Figs. 2a and 2b. The error bars and the bands 
represent solely statistical fluctuations, (a) Total projection, (b) high-energy projection, (c) low-energy projection. 



down close to the edges of the Dalitz plot [28, 29]. Furthermore, in this region, the kinematical configuration of the 
final state momenta renders the quasi two-body treatment less trustworthy as well. Finally, our model does not include 
the tensor component. Although marginal, this amplitude has a non-trivial distribution in the phase space and could 
induce sizable interference effects in our plots. In the vector channel, we find puzzling that the A'* (1410), which gives 
a sizable contribution for — > KnVr [23, 25], is hardly seen in experimental analyses of D+ — > K^K^n^. 

In conclusion, since we do not fit the Dalitz plot we think that the agreement between the model and the experimental 
data is satisfactory. 
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